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The status of the physics of heavy ion collisions is reviewed based on measurements 
over the past 6 years from the Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory. The dense nuclear matter produced in Au+Au collisions with 
nucleon-nucleon cm. energy ^sjvjv = 200 GeV at RHIC corresponds roughly to the 
density and temperature of the universe a few microseconds after the 'big-bang' and has 
been described as "a perfect liquid" of quarks and gluons, rather than the gas of free 
quarks and gluons, "the quark-gluon plasma" as originally envisaged. The measurements 
and arguments leading to this description will be presented. 

1. Introduction 

High energy nucleus- nucleus collisions provide the means of creating nuclear matter 
in conditions of extreme temperature and density CHEl. The kinetic energy of the 
incident projectiles would be dissipated in the large volume of nuclear matter in- 
volved in the reaction. The system is expected to come to equilibrium, thus heating 
and compressing the nuclear matter so that it undergoes a phase transition from 
a state of nucleons containing bound quarks and gluons to a state of deconfined 
quarks and gluons, in chemical and thermal equilibrium, covering the entire volume 
of the colliding nuclei or a volume that corresponds to many units of the charac- 
teristic length scale. This state of nuclear matter was originally given the name 
Quark Gluon Plasma(QGP) a plasma being an ionized gas. However the results 
at RHIC to be presented here indicated that instead of behaving like a gas of free 
quarks and gluons, the matter created in heavy ion collisions at nucleon-nucleon 
cm. energy ^/snn — 200 GeV appears to be more like a liquid This matter 
interacts much more strongly than originally expecte d, as elaborated in recent peer 
reviewed articles by the 4 RHIC experiments I ^ I ^ I ^ ISI^ which inspired the theorists^ 
to give it the new name "sQGP" (strongly interacting QGP). 
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Two energy regimes are discussed for the QGP At lower energies, ^/snn — 
5—10 GeV, typical of the AGS and CERN fixed target programs ^ the colliding 
nuclei are expected to stop each other, leading to a baryon-rich system. This will 
be the region of maximum baryon density. At very high energy, ^snn ?^100 — 200 
GeV, nuclei become transparent and the nuclear fragments will be well separated 
from a central region of particle production at mid-rapidity. This is the region 
of the baryon-free or gluon plasma, while in the nuclear fragmentation regions a 
baryon-rich plasma may also be formed ^1^. 

In the terminology of high energy physics, the QGP or sQGP is called a "soft" 
process, related to the QCD confinement scale 

AQj,D^(0-2GeV)-^c.lfm . (1) 

With increasing temperature, T, in analogy to increasing , the strong coupling 
constant as{T) becomes smaller, reducing the binding energy, and the string ten- 
sion, cr{T), beco mes smaller, increasing the confining radius, effectively screening 
the potentiaPl 

V{r) = -- ho-r^--— +a (2) 

3 r 6 r /I 

where ^ = IJ.{T) — l/rjj is the Debye screening mass^^. For r < l//i a quark feels 
the full color charge, but for r > 1/fJ., the quark is free of the potential, effectively 
deconfined. 

There has been considerable work over the past three decades in making quan- 



titative predictions for the QGP . The predicted transition temperature from a 
state of hadrons to the QGP varies, from Tc ~ 150 MeV at zero baryon density, 
to zero temperature at a critical baryon density roughly 1 GeV/fm^, ~ 6.5 times 
the normal density of cold nuclear matter, po = 0.14nucleons/fm^, fiB — 930 MeV, 
where fj,B is the Baryon chemical potential. A typical expected phase diagram of 
nuclear matter is shown in Fig. [T^. Not distinguished on Fig.[T^ in the hadronic 
phase are the liquid self-bound ground state of nuclear matter and the gas of free 
nucleons 

Predictions for the transition temperature for fis ^ are constrained to a 
relatively narrow range 140 < Tc < 250 MeV, while the critical energy density is 
predicted to be 5 to 20 times the normal nuclear energy density, eo = 0.14 GeV/fm'^. 
Presumably, the most accurate predictions of the pha se t ransition ar e giv en by nu- 
merical solutions of the QCD Lagrangian on a lattice see Fig. [TJjISni Here, the 
solid line indicates a first-order phase transition at larger values of /x^ > 360 ± 40 
MeV, with critical endpoint indicated by the small square, followed a smooth 
crossover for fis ^ 360 MeV. Interestingly the /is for the critical end point corre- 
sponds to y/sNN — 10 GeV at mid-rapidity, between the maximum and minimum 
nucleon-nucleon cm. energies of the AGS and CERN fixed target programs, and 
is considerably below the value of /is ~ 24 MeV for mid-rapidity at RHIC (see 
later discussion). 
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Fig. 1. a) (left) Proposed phase diagram for nuclear matter'^-': Temperature, T, vs Baryon Chemi- 
cal Potential, ^. b) (right) Lattice calculation with "2+1 staggered quarks of physical masses" 1221 



A nice feature of the search for the QGP is that it requires the integrated 
use of many disciphnes in Physics: High Energy Particle Physics, Nuclear Physics, 
Relativistic Mecha nics, Quantum Statistical Mechanics, and, recently, AdS/CFT 
string theorylI "^ l ^^ i From the point of view of an experimentalist there are two major 
questions in this field. The first is how to relate the thermodynamical properties 
(temperature, energy density, entropy, viscosity ...) of the QGP or hot nuclear 
matter to properties that can be measured in the laboratory. The second question 
is how the QGP can be detected. 

One of the major challenges in this field is to find signatures that are unique to 
the QGP so that this new state of matter can be distinguished from the "ordinary 
physics" of relativistic nuclear collisions. Another more general challenge is to find 
effects which are specific to A+A collisions, such as collective or coherent phenom- 
ena, in distinction to cases for which A+A collisions ca n be cons idered as merely 
an incoherent superposition of nucleon-nucleon collisions I^^HMini 



2. J/^f suppression — the original "gold-plated" QGP signature 

Since 1986, the 'gold-plated' signature of deconfinement was thought to be J/^* sup- 
pression. Matsui and Satz 1221 proposed that J /"^ production in A+A collisions will 
be suppressed by Debye screening of the quark color charge in the QGP. The J/"^ 
is produced when two gluons interact to produce a c, c pair which then resonates 
to form the J/\l/. In the plasma the c, c interaction is screened so that the c, c go 
their separate ways and eventually pick up other quarks at the periphery to become 
open charm. J /'^ suppression would be quite a spectacular effect since the naive 
expectation was that J/^E* production, due to the relatively large ~ 1.5 GeV scale of 
the charm quark mass, should behave like a pointlike process with production cross 
section proportional to A {BA) for p+A (B+A) minimum bias collisions, and thus 
would be enhanced relative to the total interaction cross section, which increases 
only as A^/^ . However, there were problems from the very beginning because unlike 
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Drell-Yan which exhibits an dependence in p+A coUisions, the J/^P is sup- 

pressed ^23by A*^-^^. The suppression continues as {BA)'^-^'^ for minimum bias B+A 
colhsions'241 (Fig. [2^) at the SPS until the heaviest system, Pb+Pb, where the sup- 
pression increases by ~ 25%. However the suppression is much more impressive as 
function of centrality (Fig. [2)3) . This is the CERN fixed target heavy ion program's 
main claim to fame. 
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Fig. 2. a) (left) Total cross section for J/^ production divided by AB in A+B collisions at 158— 
200A GeV EH. b) (right) Cross section for J/"!- divided by Drell-Yan as a function of centrality 
measured by Ex for 158A GeV Pb-|-Pb collisions ESI. The solid line represents the expected 
suppression for cold nuclear matter. The plot may be converted to Raa by dividing by 34. 



My summary of the different views of dilepton resonances in the High Energy 
and Relativistic Heavy loi*^ Physics communities since the mid 1980's is shown 
in Fig. H 



3. Observables in Relativistic Heavy Ion Collisions 

A schematic drawing of a collision of two relativistic Au nuclei is shown in Fig. |4l 
In the center of mass system of the nucleus-nucleus collision, the two Lorentz- 
contracted nuclei of radius R approach each other with impact parameter b. In the 
region of overlap, the "participating" nucleons interact with each other, while in 
the non-overlap region, the "spectator" nucleons simply continue on their original 
trajectories and can be measured in Zero Degree Calorimeters (ZDC), so that the 
number of participants can be determined. The degree of overlap is called the 
centrality of the collision, with b ^ being the most central and b ^ 2R, the 
most peripheral. The maximum time of overlap is tq = 2R/^c where 7 is the 
Lorentz factor and c is the velocity of light. The energy of the inelastic collision is 
predominantly dissipated by multiple particle production, where rich, the number of 
charged particles produced, is directly proportional to the number of participating 
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Success in HEP Success in RHI 
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Fig. 3. "The road to success": In High Energy Physics (left) a UAl mcasuremenl!2Si of pairs of 
muons each with px > 3 GeV/c shows two Nobel prize winning dimuon peaks and one which won 
the Wolf prize. Success for measuring these peaks in RHI physics is shown schematically on the 
right. 
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Fig. 4. a) (left) Schematic of collision of two nuclei with radius R and impact parameter b. The 
curve with the ordinate labeled dcr/dn^i^ represents the relative probability of charged particle 
multiplicity n^h which is directly proportional to the number of participating nucleons, A'^part. 



nucleons (iVpart) as sketched on Fig. H) Thus, rich in central Au+Au collisions is 
roughly A times larger than in a p-p collision, as shown in actual events from the 
STAR and PHENIX detectors at RHIC in Fig. E 

It would appear to be a daunting task to reconstruct all the particles produced 
in such events. Consequently, it is more common to use single-particle or multi- 
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Fig. 5. a) (left) A p-p collision in the STAR detector viewed along the collision axis; b) (center) 
Au+Au central collision at ^sjvjv = 200 GeV in the STAR detector; c) (right) Au+Au central 
collision at v^sjvjv = 200 GeV in the PHENIX detector. 



particle inclusive variables to analyze these reactions. For any observed particle of 
momentum p, energy E, the momentum can be resolved into transverse {pt) and 
longitudinal {pl) components; and in many cases the mass (m) of the particle can 
be determined. The longitudinal momentum is conveniently expressed in terms of 
the rapidity (y): 

.^l.f^) (3) 

\ TUT J 

coshy = E/rriT smhy = PL/friT dy = dpi^/E (4) 

where 



Tot = \/ + Pt 



and E = \Jp\+ = ^p'^ +m? (5) 

In the limit when (to <C E) the rapidity reduces to the pseudorapidity (77) 

7] = -lntane'/2 (6) 

cosh 77 = CSC smhrj = coi6 (7) 

where 9 is the polar angle of emission. The rapidity variable has the useful property 
that it is additive under a Lorentz transformation. 



4. The RHIC facility 

The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Labo- 
ratory consists of two independent superconducting acceleration/storage rings of 
3.8km circumference which cross at six interaction regions where experiments can 
be performed (see Fig.[6|). At present, RHIC can collide any nucleus from proton to 
Au with any other nucleus, preferably at equal velocity /? = w/c (7 = l/-\/l — 
with an achieved luminosity of 4 x 10^® cm~^ sec~^ for Au+Au collisions at a 
beam energy of 100 GeV/nucleon which corresponds to a nucleon-nucleon collision 
cm. energy ^ snn = 200 GeV. For p-p collisions, the beams are always polarized, 
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Polarized Proton Collisions at BNL 




Fig. 6. Schematic of RHIC accelerator with emphasis on equipment for polarized proton acceler- 
ation and storage. For Heavy Ions, injection to the booster is via transfer line from the Tandem 
Van de Graaf. 

with a design polarization of 70% per beam, with 65% achieved at present with 
a maximum luminosity of 2 x 10^^ cm'^ sec-^ ai y/s = 200 GeV. The maximum 
cm. energy for p-p collisions is -^/s — 500 GeV, with a luminosity of > 10^^ cm~^ 
sec~^. In general, the luminosity in AB collisions is proportional to ^"^ /{AB) from 
the injection energy of ^ 24 Z/A GeV to the maximum energy of 250 Z/A GeV 
and is proportional to 7'^ — 7^ for proposed deceleration below the injection energy 
to allow Au+Au collisions with ^snn as low as > 5 GeV to access the region of 
the possible critical point (Fig. [Ha). In the future, electron cooling will increase the 
luminosity by a factor > 20 and a new Electron Beam Ion Source (EBIS) to replace 
the Tandem Van de Graaf and transfer line will allow nuclei up to Uranium A 
summary of the different species, energies and data collected in the first six years 
of RHIC operation is given in Table [TJ 

4.1. The Experiments 

There are two major experiments at the RHIC heavy ion program, STAR ^ and 
PHENIxl^, and two smaller experiments. The two small experiments, which have 
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Table 1. Summary of RHIC runs: Date, conditions, Integrated Luminosity, number of scanned 
interactions (A^tot int.), equivalent p-p Luminosity for Hard-Scattering, Data Archived. 



Run 


Year 


Species 


y/SNN 


J Ldt 




Ntot 


p-p Equiv. 


Data 








GeV 






int. 


ABx J Ldt 


Size 


01 


2000 


Au+Au 


130 


1 fih- 


-1 


lOM 


0.04 pb-i 


3 TB 


02 


2001-2 


Au+Au 


200 


24 ^b- 


-1 


170M 


1.0 pb-i 


10 TB 






P+P 


200 


0.15 pb- 


-1 


3.7G 


0.15 pb-i 


20 TB 


03 


2002-3 


d+Au 


200 


2.74 nb" 


-1 


5.5G 


1.1 pb-i 


46 TB 






P+P 


200 


0.35 pb- 


-1 


6.6G 


0.35 pb-i 


35 TB 


04 


2003-4 


Au+Au 


200 


241 ^b- 


-1 


1.5G 


10.0 pb-i 


270 TB 






Au+Au 


62 


9 fih' 


-1 


58M 


0.36 pb-i 


10 TB 


05 


2004-5 


Cu+Cu 


200 


3 nb- 


-1 


8.6G 


11.9 pb-i 


173 TB 






Cu+Cu 


62 


0.19 nb" 


-1 


0.6G 


0.8 pb-i 


48 TB 






Cu+Cu 


22.5 


2.7 Aib- 


-1 


9M 


0.01 pb-i 


1 TB 






P+P 


200 


3.8 pb- 


'1 


85B 


3.8 pb-i 


270 TB 


06 


2006 


P+P 


200 


10.7 pb- 


-1 


230B 


10.7 pb-i 


310 TB 






P+P 


62 


0.1 pb- 


-1 


28B 


0.1 pb-i 


25 TB 



now completed their programs, are PHOBOS'SJ, which emphasizes charged particle 
detection over the full phase space using Si detectors, and BRAHMS ^, with two 
small-aperture, high precision moveable spectrometers, one at mid-rapidity and one 
capable of moving as far forward as 2° {\ri\ > 4). 

STAR, which emphasizes hadron physics, is most like a conventional hadron- 
coUider detector, a TPC covering the full azimuth over ±1 unit of pseudorapidity, for 
the purpose of charged particle tracking in a magnetic field of 0.5 Tesla. The TPC is 
surrounded by a system of Time of Flight counters, for particle identification, and 
a moderate resolution {15%/^/E) electromagnetic calorimeter, for measuring tt" 
production and charged-neutral energy correlations. The detector is completed by 
a Silicon Drift Vertex Tracker, for measurements of Hyperons, and possible TPC's 
external to the magnet, for tracking at small angles 2.0 < \ri\ < 4.5. 

PHENIX, a very high granularity, high resolution detector for leptons and 
photons emerging from the Quark Gluon Plasma (QGP), emphasizes the abil- 
ity to run at the highest luminosities with very selective triggers to find rare 
events, particularly J/'^, open charm and direct photon production. PHENIX has 
a highly instrumented electron, photon and charged hadron spectrometer, in the 
central region |?7| < 0.35, with full azimuth di-muon measurement in two endcaps, 
1.15 < |?7| < 2.35. The electron/photon central spectrometer emphasizes electron 
identification at the trigger level, with RICH, TRD and EM calorimetry. The EM 
calorimeter, with energy resolution as/E = ^%/^/E{GeV), also serves as an ex- 
cellent photon and 7r° trigger because of its 5 by 5 cm segmentation at 5.1 m. 
The central spectrometer consists of two arms, each subtending 90° in azimuth (<I>) 
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Fig. 7. a (left) Schematic of 7r± , e± and 7 in PHENIX, with ElectroMagnetic Calorimter (EMCal) 
and Ring Imaging Cherenkov Counter (RICH), b) (right) Energy/momentum for all charged 
particles detected in the EMcal with and without a RICH signal. 

and ±0.35 units in pseudorapidity (77). The total coverage is 1/2 of the azimuth — 
however, the two arms are not back-to-back: the gap between the edges of the two 
90° arms is 67.5° on one side and hence 112.5° on the other. The charged parti- 
cle momentum resolution is 1% at 5 GeV/c, and charged hadron identification is 
provided by TOF(lOOps) for 1/3 of the azimuth of one arm. 

Rather than show a diagram of PHENIX, which is a very non-conventional 
collider detector [2Si^ the principles of the design of PHENIX will be illustrated 
(see Fig. [71). In order to detect electrons in hadron collisions, where the typical 
ratio of b^/tt^ is known to be of order of 10^^ for prompt electrons (from charm) 
at 23.5 < y/s < 62.4 GeV (ISR energies), one must plan on a charged pion 
rejection of > 10^. PHENIX decided to use a Ring Imaging Cerenkov counter 
(RICH) in combination with an Electromagnetic Calorimeter (EMCal) to achieve 
this rejection. There is also the issue of the huge background of from internal 
and external conversion of the photons from 7 -I- 7 decay or from direct photon 

production, which must be measured and understood to high precision. The axial 
magnetic field in PHENIX is designed with the possibility of zero magnetic field on 
the axis — so that the lower energy member of an asymmetric e~^e~ conversion pair 
does not curl up and get lost — and with a minimum of material in the aperture 
(0.4% of a radiation length) to avoid external conversions. The EMCal is crucial for 
electron identification and hadron rejection. The EMCal measures the energy of 7 
and and reconstructs tt" from 2 photons. The high granularity of the individual 
EMCal towers is {drj x 6(j)) ^ (0.01 x 0.01), which allows the two photons from a 
7r° to be resolved from a single photon cluster for values of pr up to > 30 GeV/c. 
It measures a decent time of flight (TOE), 0.3 nanoseconds over 5 meters, allowing 
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photon and charged particle identification. The high precision TOF over part of 
the aperture allows improved charged hadron identification. Electrons are identified 
by a count in the RICH and matching Energy and momentum (E/p), where the 
momentum is measured by track chambers in a magnetic field. Charged hadrons 
deposit only minimum ionization in the EMCal (~ 0.3 GeV), or higher if they 
interact, and don't count in the RICH (tt^ threshold 4.7 GeV/c). Thus, requiring 
a RICH signal rejects all charged hadrons with p < 4.7 GeV/c, leaving only 
as indicated by the E/p = 1 peak in Fig. [7]-right. tt^ can be identified above 4.7 
GeV/c by a RICH signal together with an EMCal hit with minimum ionization or 
greater. It is amusing to realize that once you decide to measure electrons, you must 
make an excellent tt" measurement to understand the background, and this implies 
a detector which can measure and identify almost all particles, as exemplified by 
the PHENIX central spectrometer. 



5. Inclusive and semi-inclusive particle production 

A single particle "inclusive" reaction involves the measurement of just one particle 
coming out of a reaction, 

a + 6 — > c + anything . 

The terminology 1^ comes from the fact that all final states with the particle c are 
summed over, or included. A "semi-inclusive" reactiorP^ refers to the measurement 
of all events of a given topology or class, e.g. 

a + — > ni particles in class 1 + anything , 

where "centrality" is the most common class in relativistic heavy ion collisions. 

Measurements are presented in terms of the (Lorentz) invariant single particle 
inclusive differential cross section (or Yield per event in the class if semi- inclusive): 

Ed^^ = ^ 1. f (p^ y) (10) 

dp^ PTdprdydcj) 2n ' ' 
where y is the rapidity, px is the transverse momentum, and (j) is the azimuth of 
the particle (see Fig. |8l). The average transverse momentum, (pr), or the mean 
transverse kinetic energy, (mr) — m, or the asymptotic slope are taken as measures 
of the temperature, T, of the reaction. 

It is important to be aware that the integral of the single particle inclusive cross 
section over all the variables is not equal to cr/ the interaction cross section, but 
rather is equal to the mean multiplicity times the interaction cross section: (n) x aj. 
Hence the mean multiplicity per interaction is 



(n) = — / ^dydpT Pr f (pt, v) = — f = [ 

(71 J 27r cri J dy J 



dyp{y) , (11) 



where the terminology for the multiplicity density in rapidity is {l/ai)da/dy — 
p{y) = dn/dy for identified particles (m known), dn/drj for non-identified particles 
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Fig. 8. a) (left) Se mi-in clusive invariant pT spectra for vr^, K^, in Au+Au collisions at 
y/sjVJV = 200 GeV IS132] b) (right) (pt) of positive particles as a function of centrality (A^part) 
from the same data. 



(m unknown, assumed massless) . The total charged particle multiplicity is taken as 
a measure of the total entropy, s and dn/drj is taken as a measure of the entropy 
density in restricted intervals of rapidity. 



5.1. The rapidity density and the rapidity plateau. 

At RHIC, the reference frame of the detectors is designed to be the nucleon-nucleon 
cm. frame in which the two nuclei approach each other with the same 7. The 
nucleon-nucleon cm. energy is denoted ^/snn, and the total cm. energy is ^/s = 
^ ■ \/snn for symmetric A-|- A collisions. The colliding nucleons approach each other 
with energy y^ijvjv/2 and equal and opposite momenta. The rapidity of the nucleon- 
nucleon center of mass (mid-rapidity) is y*^™ — unn — 0, and the projectile and 
target nucleons are at equal and opposite rapidities: 

yP™j = _ytarg<=t = cOsh-1 = ybca„Y2 ^ (12) 

where mjv =931 MeV is the mass of nucleon bound in a nucleus. 

The shape and evolution with of the charged particle density in rapidity, 
dn/dy, in p-p, p-l-A and A-fA collisions all follow a similar trend (Fig [9^) and 
provide a graphic description of high energy collisions. Regions of nuclear fragmen- 
tation take up the first 1-2 units around the projectile and target rapidity and if the 
center-of-mass energy is sufficiently high, a central plateau is exhibited. The dis- 
tributions increase in width with increasing ^^snn but by a smaller amount than 

''The Laboratory reference frame used in fixed target experiments is shifted by y''°^™/2 from the 
cm. system such that yf'<^°^ = 0, yf"' = y^<=^'-^\ y'j^^^ = y'^^^-^/2. 
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Fig. 9. a) (top right) dn/drj for A+A and d+Au collisions at RHIC as a function of y'ijvjv'^^. b) 
(bottom left) dn/dy\p - dn/dy\p at AGS, SPS and RHIC, yP^j = 1.6, 2.9, 5.4Eil. 



the increase in y^°^™ and show a smah decrease in width with increasing centrahty. 
dn/drj increases with increasing centrahty, ^/sNN and A in A+A coUisions. In the 
asymmetric d+Au cohision, dn/drj in the target rapidity region is larger than in 
the projectile region, but not by much, only about 50%. Also the nuclear trans- 
parency is evident, there is no reduction of particles at the projectile rapidity with 
increasing centrahty. 

Subtleties of the dist ribu tions in A+A collisions become apparent when iden- 
tified particles are used In Fig. [91d, the difference of dn/dy for protons and 
anti-protons, i.e net-protons is shown as a function of cm. energy, y^SNN = 5 
(AGS, Au+Au), 17 (SPS, Pb+Pb), 200 (RHIC, Au+Au) GeV, yP^j =1.6, 2.9, 5.4. 
As ^/s'nn is reduced, stopping of the participating nucleons is indicated by the 
nuclear fragmentation peak moving from the fragmentation region (not visible for 
RHIC) to mid-rapidity. 
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Fig. 10. a) (left) PHENIxE3_Bt distribution in Au+Au collisions at v^sjvjv = 200 GeV showing 
centrality selections of 0-5%, 5-10%, 10- 15%. 15-20%. b) (right) Event-by-event average By/cluster 
in the EMCal for 20-30% centrality ES] (data points), random baseline (dotted line). 



5.2. Centrality Measurement, Fluctuations, Bjorken Energy 
Density 

Another variable, closely related to multiplicity, is the transverse energy density in 
rapidity or dEx/dy ~ (pr) x dn/dy, usually measured in calorimeters by summing 
over all particles on an event in a fixed but relatively large solid angle Et = 
J2i EiSindi. dEx/dy is thought to be related to the co- moving energy density in 
a longitudinal expansion IS^IHl^ ^nd taken by all experiments as a measure of the 
energy density in space e: 

where Tp, the formation time, is usually taken as 1 fm/c, irR^ is the effective area 
of the collision, and d (Et) / dy is the co-moving energy density. 

Besides the Bjorken energy density, the importance of Et distributions in RHl 
collisions (see Fig. llOb ) is that they are sensitive primarily to the nuclear geometry 
of the reaction, and hence can be used to measure the centrality of individual inter- 
actions on an event-by-event basis. However, in PHENIX, centrality is determined 
far from mid-rapidity to avoid possibly biasing any of other the mid-rapidity mea- 
surements. The centrality is determined by the upper percentiles of the charged 
multiplicity distribution in Beam Beam Counters (BBC) which cover the region 
3.0 < \r}\ < 3.9 in coincidence with Zero Degree Calorimeters located at |?7| > 6. 
The shift of the Et distributions to higher values with increasing centrality (smaller 
upper-percentiles) is nicely illustrated in Fig. [TOk . For 0-5% centrality the Bjorken 
energy density is 5.5 GeV/fm^, well above the value of 1 GeV/fm'^ nominally envi- 
sioned for QGP formation. 

Fluctuations can also be studied. For instance, is the shape of the upper edge of 
the centrality selected Et distributions in Fig. llOb random, or is it evidence of non- 
random fluctuations? The event-by-event distribution Mej, of the 
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cluster, which is closely related to Et and also to Mp^ , the event-by-event average 
Pt of charged particles, is shown in Fig. llOb ^^Sl. The non-random fluctuation is 
measured by comparison to a random baseline from mixed events. Mp^ is defined: 

1 " 1 

Mp^^W=-Y.PT.'^ -Etc , (14) 

i=l 

where Etc is the analog of Et for charged particles. PHENIX has shown that 
the non-random Mp^ fluctuations are of the order of 1%. This places severely small 
limits on the critical fluctuations expected for a flrst order phase transition but is 
consistent with the lattice QCD prediction of a smooth crossover at RHIC energies 
(Fig. Ed). 



5.3. Particle abundances, thermal/ chemical equilibrium 

As shown in Fig. [SId, the (pt) of tt^, K^, p, p, at RHIC, increases smoothly from 
peripheral to central Au-|-Au collisions, and as in pp collisions increases with in- 
creasing mass as would be expected for a thermal distribution (Eq. [15]) : 

= ^''^ = ^ ^ e-^/^ . (15) 
dpLPxdpT dpLmTdrriT e^l^ ± 1 

Since E = rriT coshy, a signal of thermal production is that the pt and mass depen- 
dence of the cross section are not independent but depend only on the transverse 
mass rriT. Thus, the inverse-slope of the rriT distribution at mid-rapidity should 
represent the temperature. However, the mT-scaling property of thermal distribu- 
tions is not exactly exhibited by the data (see Fig. [Tib ) . The inverse slopes increase 
linearly with rest mass, by a larger amount with increasing centrality, which is ev- 
idence for collecti ve motio n ('radial flow') and is seen in RHl collisions at AGsl^U^ 
SPS^^and RHIC IS[ 6 | 7 | 8 | energies. Particles (or partons) which travel with a trans- 
verse flow velocity Pt acquire kinetic energy in addition to the thermal energy so 
that the inverse slope should increase linearly with the rest mass, T —f Tq + jt wq, 
as illustrated by the lines on Fig. [TTk . The effect is primarily at low pT where the 
slope flattens with increasing centrality as illustrated in Fig. [TTb . 

The semi-inclusive ratios of different particle abundances also vary smoothly 
as a function of centrality in Au+Au collisions at RHIC (see Fig. with a 

considerably larger increase in production than p, p production relative to tt^ . 
However the p/n~^ and p/tt" ratios as a function of pt (Fig. 112b ) show a dramatic 
increase as a function of centrality at RHIC which was totally unexpected and 
is still not fully understood (see below). The ratios of particle abundances (which 
are dominated by low pT particles) for central Au-|-Au collisions at RHIC, even 
for strange and multi-strange particles, are well described (Fig. [T5b ) by flts to a 
thermal distribution, 
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Fig. 11. a) (left) Inverse slope of the my distribution l^^^l (colloquially called T) for vr^, , p, 
p in Au+Au collisions for three centralities at .y/sjvjv = 200 GeV. b) (right) Inviarant yield of p 
and p as a function of centrality scaled (by the number of binary collisions N^^n) to lie on top of 
each other for pT>2 GeV/cHSl 
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Fig. 12. a) (left) K^/tt^ and p^/n^ ratios as a function of centrality (A^part) in_^/sjvjv Au+Au 
coUisionsE! b) (right) p/ir and p/ir ratio as a function of px from the same 
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points arc for vr^ ('""'' )i respectively. 



with similar expressions for strange particles, where fj,B and fis ^re chemical poten- 
tials associated with each conserved quantity: baryon number (/is) and strangeness 
(us)- This should not be very surprising as the particle abundances in A+A colli- 
sions at SPS and AGS energies and in p-pS5!and e~^e~ collisions'^ are also well 
described by the same thermal model to such an extent that the chemical freezeout 
temperature Tch as a function of /is (from Eq. I16p could be derived, which looks 
suspiciously like a phase diagram (Fig. [T3b). Of course, as the thermal equilibrium 
properties of the QGP are the subject of interest, it is important to understand how 
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Fig. 13. a) (left) Ratios of pr-integrated mid-rapidity yields for different hadron species measured 
in Au+Au central central collisions at yTf^ = 200 GeV 43, BRAHMS(circles) and PHOBOS 
(crosses) data are shown as well as the PHENIX (triangles) and S TAR (stars) data indicated, b) 



(right) T|,jj versus for thermal model fits as a function of y^Tl 



EH 



where the line represents a 



parameterization of the freeze-out curve, the temperature where chemical equilibrium is achieved 
so that the particle abundances are fixed. 



or if the thermal properties of the observed hadrons relate to the thermal proper- 
ties of quarks and gluons in the QGP before getting too excited. Also, while these 
thermal models appear to be simple they have many technical details which are 
beyond the scope of this article. 



6. Flov^r 

A distinguishing feature of A-t-A collisions compared to either p-p or p-l-A collisions 
is the collective flow observed. This effect is seen over the full range of energies 
studied in heavy ion collisions, from incident kinetic energy of lOOA MeV to cm. 
energy of ^snn ~ 200 GeVl^. Collective flow, or simply flow, is a collective effect 
which can not be obtained from a superposition of independent N-N collisions. 

Immediately after an A-l-A collision, the overlap region defined by the nuclear 
geometry is almond shaped (see Fig [Til) with the shortest axis along the impact pa- 
rameter vector. Due to the reaction plane breaking the 4> symmetry of the problem, 
the semi-incl usiv e single particle spectrum is modified from Eq.fTUlby an expansion 
in harmonics o f the azimuthal angle of the particle with respect to the reaction 
plane, (t> — where the angle of the reaction plane '^r is defined to be along 
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Plane \^ 






Fig. 14. a) (left) Almond shaped overlap zone generated just after an A+A collision where the 
incident nuclei are moving along the ±2 axis, and the reaction plane, which by definition cont ains 
the impact parameter vector (along the x axis) Thanks to Masashi Kaneta for the figure HHl. 
b) (right) View of the collision down the z axis: (top) spatial distribution (bottom) momentum 
distribution after elliptic flow {v2) develops 



the impact parameter vector, the x axis in Fig. 1141 



dp^ 



prdpTdydcf) 
d^N 



(17) 



-[1 + 2vi COs{(j> — + COS 2(0 — $/ 



2^ pTdpTdy 

The expansion parameter vi is cahed the directed flow and V2 the eUiptical flow. If 
no collective behavior takes place, i.e. the interaction is merely a superposition of in- 
dependent nucleon-nucleon collisions, then the outgoing momentum distribution of 
the particles would be isotropic in azimuth. However, since the leading participating 
nucleons in the forward region +z (Fig. [Mb ) will interact with many other nucleons 
in the "almond", they will be pushed away from the rest of the participants, into 
the -l-a; direction, while the —z going participants are pushed towards — x. This is 
what causes the directed flow, vi, which was discovered at the Bevalac^^and is 
clearly sensitive to the Equation of State. For instance if one imagines the almond 
to be composed of billiard balls requiring lots of pressure for a small deformation 
(hard EOS) a larger vi would result than if the almond suddenly melts, perhaps 
turning into a 'perfect fluid', with a much softer EOS^^. 

The same principles apply to W2, the parameter of cos 2(0 — ^r), which (un- 
like vi) doesn't change sign with rapidity, and hence is non-zero at midrapidity. If 
thermal equilibrium is reached, then the pressure gradient is directed mainly along 
the direction of the impact parameter {x axis in Fig. 114b ) and collective flow de- 
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velops along this direction. II all the particles are approximately at rest in the fluid 
and thus move with the fluid velocity, the transverse momentum distribution will 
reflect the fluid profile. Hence the anisotropic spatial distribution is c arri ed over to 
an anisotropic momentum distribution through the pressure gradient 

It is important to emphasize that the spatial anisotropy turns into an mo- 
mentum anisotropy only if the outgoing particles (or partons) interact with each 
other Thus the momentum anisotropy is proportional to the spatial anisotropy 
of the almond, represented by the eccentricity, e = {Ry — R'^)/{Ry + R^) ~ 
{Ry — Rx)/{Ry + Rx), at the time (io) of thermalization. This is due to the fact 
that the meannumber of scatterings in the transverse plane is different along the 
X and y axes I 54 | 55 | 50 | rpj^^ mean number of scatterings is proportional to the par- 
ticle density, p = (l/nRxRy) dn/dy (similar to Eg. I13p times the interaction cross 
section {a) times the distance traversed: 

1 dn 1 dn I dn 

V2(xRy<7 Rxa — — (X £g — — , (18) 

irRxRy dy irRxRy dy nRxRy dy 



where Rx = a/ (x^), Ry = ^/{y^■ 




"ch/"max p, (GeV/c) 



Fig. 15. a) (left) V2 a s fu nction of centrality (measured as a fraction of the total number of charged 
particles nch/^-max) The boxes represent the expected hydro-limit with V2/e = 0.19 (lower 
edge) and 0.25 (upper edge), b) (right) V2 as a function of for identified particles in minimum 
bias Au-|-Au collisions at ^smn = 200 GeV together with a hydro calculation EH. 



Since the eccentricity s is much larger for peripheral than for central colli- 
sions, the dep endence of V2 on centrality should exhibit a characteristic shape (see 
Fig. [l3b ) . This was one of the first publications from RHIC and showed that 
V2 was surprisingly large and near the hydro-limits. The hydro-limits indicated are 
for full thermalization of the system at the value of e given by the initial nuclear 
geometry of the almond at the time of overlap. If the system doesn't thermalize 
rapidly, the fiow tends to vanish because the eccentricity reduces as the system 
expands ° 

Another surprise!^ (Fig. ITSb) was that the V2 followed the hydro prediction 
out to ~ 2 GeV/c and then plateaued at a constant value to much higher pT- 
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This was one of the principal arguments for the "perfect fluid" because any modest 
value of viscosity would cause the V2 to decrease towards zero near pT ~ 1.7 
GeV/c (see Fig.^, below). 
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Fig. 16. a) (left) 1)2/1 vs pt/tl for identified particles, where n is the number of constituent quarks, 
together with a polynomial fitEZl. Bottom panel shows devi ations from the fit. b) (right) V2/n vs 
transverse kinetic energy per constituent quark KEx/n^^ 



As hydrodynamics appears to work in both p-p and A+A collisions, and col- 
lective flow is observed in A+A collision over the full range of energies studied, a 
key question is what is flowing at RHIC and is it qualitatively different from the 
flow observed at lower -y/sjvjv ? One interesting proposal in this regard is that the 
constituent quarks flow , so that the flow should be proportional to the number of 
constituent quarks Uq, in which case vijuq as a function of px/nq would represent 
the constituent quark flow as a function of constituen t qu ark transverse momentum 
and would be universal. Interestingly, the RHIC data^^ (Fig.flBk) seem to support 
this picture, although the fact that the 7r+ + 7r~ deviate most from the universal 
curve should raise some suspicions as the pion is the only particle whose mass is 
much less than that of its constituent quarks. However, in relativistic hydrodynam- 
ics, at mid-rapidity, the transverse kinetic energy, — toq = (7t ~ 1)^0 = KEt, 
rather than px is the relevant variable, an d in fact vijuq as a function of KEt jriq 
seems to exhibit nearly perfect scaling 1^ (Fig. [16b). 

Another striking hint as to what is flowing at RHIC is given i n Fig. [T7b where 
charm particles dete cted by their large semi-leptonic decay^^^^^ exhibit the same 
vi as other particles l ^^ l ^^ l pgj. ^j-^g charm particles, the vi of the decay electrons 
follows the V2 of the D(cd^ cu) mesons but due to their different masses the c 

and u, d quarks have different momenta for the same velocity required for formation 
by coalescence, so the V2 of the D mesons is reduced to that of the light quark at 



February 8, 2008 21:56 WSPC/INSTRUCTION FILE MJT-BufFalo-ws- 
ijmpe 



20 M. J. TANNENBAUM 




0.2 0.4 0.6 0.8 



1.2 1.4 1.6 



0.15 
0.1 
0.05 


-0.05 



: Au+Au@\|iJjij=200GeV 


without charm flow 


~Min.Bias 






- non photonic e*+e" 


with charm flow 


— PHENIX preliminary 






r rf^ * 


f 1 






. . . 1 . . . . 1 . . . . 1 . . . . 1 , , , . 1 . . . . 1 , , , , 1 . . . . 1 


, , , ,\ , , , , 



Fig. 17. a) (loft) V2 as a function of pj^ in a hydro calculation ifi^J for mid-central collisions for 
different values of Fs/ro, the "sound attenuation length" which is zero for a "perfect fluid" and 
increases linearly with the viscosity, b) (right) V2 of non-photonic e+ -|- e~ from semi-leptonic 
heavy flavor decay ED. The solid cu rve is if both charm and light quarks flow, while the dashed 
curve is if only the light quark flows ESI. 



lower pT if the c quark itself does not flow. The data (Fig. flTb) favor the flow of 
the c quark, but clearly lots more work remains to be done to improve both the 
measurement and the theory. However since the 4>{ss) meson also flows at RHIcI^ 
and follows the constituent quark scaling rule, this clearly indicates that the flow 
at RHIC is partonic because the hadronic interaction cross section of the meson 
is much smaller than for other hadrons 70 71 72 _ 



7. Jet Quenching 

A new tool for probing the color response function of the medium was developed 
in the early 1990's^^and given a flrm basis in QCD^^just before RHIC turned 
on. In the initial collision of the two nuclei in an A-l-A collision, when the Lorentz 
contracted nuclei are overlapped, hard scatterings can occur which produce pairs 
of outgoing high pT partons with their color charge fully exposed as they travel 
through the medium before they fragment to jets of particles. If the medium also 
has a large density of color charges, then the partons are predicted to lose energy 
by 'coherent' (LPM) gluon bremsstrahlung which is sensitive to the properties of 
the medium. This leads to a reduction in the pT of the partons and their fragments 
and hence a reduction in the number of partons or fragments at a given pT, which 
is called jet quenching. The effect should be absent in p-l-A or d-f A collisions due 
to the lack of a medium produced. 

One of the r naj or , arguably the major discovery at RHIC, was the observation 
of jet quenching ''^^^1^ by the suppression of tt" and non-identified charged hadrons 
in Au-f Au collisions at mid-rapidity for large transverse momenta, > 2 GeV/c. 
In p-p collisions, particles with pt > 2 GeV/c at mid-rapidity (perpendicular to 
the collision axis) are produced from states with two roughly back-to-back jets 
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which are the result of hard-scattering of the constituents of the nucleon (current- 
quarks and gluons) as described by QCD. The su ppression of high pT particles 
was observed by all four RHIC experiments I ^ I 5 | 6 | 7 | jg ^^yi calibrated by using 
measurements in p-p, d-|-Au and Au-|-Au collisions at ^/snn = 200 GeV. The px 
spectra of tt'^ from p-p collisions at ^/s = 200 GeVlSUand Au-|-Au central collisions 
at y^SNN — 200 GeV are shown in Fig. [TH Note that the p-p measurements 




Fig. 18. a) (left) Invariant cross section EU of tt*^ at mid-rapid ity from p-p collisions at -^/ijvjv = 
200 GeV together with NLO pQCD predictions from VogelsanglS2lwith two different fragmentation 
functions, b) (right) p-p data from a) multiplied by (Ta a ) f or Au+Au central collisions (1-10%) 
plotted on a log-log scale together with the measured l^il semi-inclusive tt" invariant yield in 
Au+Au central (0-10%) collisions at yijvjv = 200 GeV. 



for pt > S GeV/c are in excellent agreement with a pQCD calculation and that 
the pt spectra in both p-p and central Au+Au follow a pure power law for px > 3 
GeV/c. 

Since hard scattering is point- like, with distance scale 1/pt < 0.1 fm, the cross 
section in p+A (B+A) collisions, compared to p-p, should be simply proportional 
to the relative number of possible point-like encounters a factor of A {BA) 
for p-|-A (B-l-A) minimum bias collisions. For semi- inclusive reactions in centrality 
class / at impact parameter b, the scaling is proportional to (Tab) f, the overlap 
integral of the nuclear thickness functions averaged over the centrality class 
/. Since {Tab) f = {Ncou) f /(^nn , where (Ncou) f is the average number of binary 
nucleon- nucleon inelastic collisions with cross section aNN for the centrality class 
/, point-like scaling is often called binary-collision (or Afco;/)-scaling. However, the 
scaling has nothing to do with the inelastic hadronic collision probability, it is 
proportional only to the geometrical factor {Tab) f- 

Effects of the nuclear medium, either in the initial or final state, may modify the 
point-like scaling. This is shown rather dramatically in Fig. 118b . where the Au+Au 
central (0-10%) collision data are suppressed relative to the scaled p-p data by a 



February 8, 2008 21:56 WSPC/INSTRUCTION FILE MJT-BufFalo-ws- 
ijmpe 



22 M. J. TANNENBAUM 



factor of ~ 4 — 5 for > 3 GeV/c. While the suppression of tt'^ at a given pT in 
Au+Au compared to the scaled p-p spectrum may be imagined as a loss of these 
particles due to, for instance, the stopping or absorption of a certain fraction of the 
parent partons in an opaque medium, it is evident from Fig. 118b that an equally 
valid representation can be given by a downshift of the scaled p-p spectrum due to, 
for instance, the energy loss of the parent partons in the medium. 

A quantitative evaluation of the suppression is made using the "nuclear modifi- 
cation factor" , Rab, the ratio of the measured semi-inclusive yield to the point-like 
scaled p-p cross section: 



Rab — 



dN^B 



dN^B 



{Tab)} X dcr^^ {Ncon)f x dN^^ 



(19) 



where dN^^ is the differential yield of a point-like process P in an ^ -I- B collision 
and da^j^ is 1 
scaling, Rab 



and da^j^ is the cross section of P in an NN (usually p-p) collision. For point-like 
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Fig. 19. a) (left) N ucle ar modification factor Raa{pt) for tt" in central and peripheral Au+Au at 
y'sjvjv = 200 GeVl53!. The error bars include all point-to-point experimental (p-p, Au+Au) errors. 
The shaded bands represent the fractional systematic uncertainties for each centrality which can 
move all the points at that centrality up and d own t ogether, b) (right) Compilation of Raa^Pt) 
for tt" production in A+A collisions from Refs. 



Fig.flQk shows that there is no suppression of vr" in Au-I-Au peripheral collisions. 
However, the suppression in Au-f Au central collisions, although quite impressive in 
its own right, is even more dramatic when compared to previous data. All previous 
measurements of nuclear effects at high > 2 GeV/c in p-f A and A+A collisions 
at lower ^Jsnn have given results which are larger than point-like scaling (Fig. 119b ). 
a situation called the 'Cronin Effect' and thought to be due to the multiple scat- 
tering of the incident partons in the nuclear matter before the hard-collision i 89 | 90l _ 
The suppression observed at RHIC is a totally new effect. 

Naturally, the first question asked about the RHIC suppression was whether 
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it is an initial state effect, produced, for instance, by 'shadowing' of the structure 
functions in nuclei, or a final state effect produced by the medium. Although orig- 
inally answered by all 4 R HIC experiments by the observation of no suppression 
in d+Au collisions I 5 | 6 | 7 | 8 |^ clearer answer came from later measurements of 
QCD hard-photon production {gq jq) i^Sj^ Fi g. [20k. and the total yield of charm 
particles {gg cc) deduced from measurements l "*^ ! "^ ! of non-photonic (Fig. l2ip 
in p-p and Au+Au collisions. Both these reactions are sensitive to the same initial 



PHENIX Au+Au (central collisions): 




r I nl , , , ^ , , , , I , , .". , , , , I 

2 4 6 8 10 12 14 16 18 20 5 10 15 20 

p,(GeV/c) "^^'"'"^ 



Fig. 20. a) (left) Nuclear modifi cati on factor, Raa for direct photons, and rj in Au-|-Au central 
collisions at y^sjviv = 200 GeV R^a for identified tt'' and-non identified charged hadrons 

{h+ +h-)/2 for central (0-10%) Au-f Au coUisions at ^t^Wn = 200 GeV EH. 



state partons as 7r° production, but the photons, which are themselves elementary 
constituents which participate in and emerge directly from the hard-scattering, do 
not interact with the final state medium; and the number of cc pairs produced 
does not depend on whether or not the c or c later interact with the medium. 
Hence the fact that Raa = 1 as a function of centrality for these two reactions, 
which is dramatically different from the suppression of tt*^, indicates that the tt*^ 
suppression is produced by the interaction of the outgoing hard-scattered parton 
with the medium, losing energy before it fragments into a tt*^. This conclusion is 
reinforced by the equal suppression of ?7-mesons and tt*^ which are both fragments 
of jets. The curve on Fig. [20k shows a theoretical prediction from a model 

of parton energy loss. The model assumes an inital parton density dN/dy = 1200, 
which corresponds to an energy density of approximately 15 GeV/fm"^. The theory 
curve appears to show a reduction in suppression with increasing pT while the data 
appear to be flat to within the errors which clearly could still be improved. 

A major surprise was revealed by the measurement of the pt dependence of 
non-photonic from charm particles, which exhibits a suppression comparable 
to tt" (Fig. 121b ) , indicating a surprisingly strong interaction with the medium 
This is is a very recent exciting result which is causing a detailed reevaluation of 
the theoretical models because the heavy c quark (and the heavier h quark) were 
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Fig. 21. a) (left) Invariant cross section 1^21 for prompt (non-photonic) electrons in p-p collisions 
at ^/s = 200 GeV compared to a theoretical calculation (see Ref. I63t which shows contributions 
from both c and b quarks, which depend on py. b) (right) Raa a-s a function of centrality (A^part) 
for the total yield of from charm (py > 0.3) GeV/c, compared to the suppression of t he 
yield at large pT > 3.0 GeV/c which is comparable to that of tt" with {pT > 4 GeV/c)ISll. 



predicted to lose much less energy to the medium than light quarks and gluons 

Another important result, as indicated in Fig. [^Db by the difference of Raa of 
unidentified charged hadrons + h^) and 7r° in the range 1.5 < pt < 5 GeV/c, 
is that protons and anti-protons are not suppressed at all in this pt range, but 
follow point-like scaling unlike the tt^ which are suppressed, which results in an 
anomalously \a,vgep^/-n^ ratio (Fig. [T2b). Thep^/ir^ ratio returns to the normal 
value consistent with jet-fragmentation forp^ ^ 5 GeV/c. (The Ncoii scaling of the 
Pt spectra can be seen in Fig. [TTb). This observation, called the 'baryon anomaly', 
remains unexplained. 



7.1. Jet properties from two-particle correlations 

The study of jet properties via two-particle correlations, pioneered at the CERN- 
ISR^SHI^ is used at RHIC rather than reconstructing jets from all particles within 
a cone of size Ar = •\/(A^p + (A^p). This is because the large multiplicity in 
A-l-A collisions (recall Fig. [lb, c) results in a huge energy 7r(Ar)^ x 5^^^ ~ 375 
GeV in the standard cone, Ar — 1, for Au-fAu central collisions. Elliptical flow 
further complicates the jet measurement because the large non-jet background is 
modulated by cos 20, which is comparable in width to the jets studied so far. In 
Fig. [22^, the conditional probability of finding an associated charged particle with 
2 < PTa ^ PTt per trigger charged particle with 6 < pTt < 8 GeV/c is shown for 
p-p and Au+Au central collisions, where the p-p data have been added to the large 



February 8, 2008 21:56 WSPC/INSTRUCTION FILE MJT-BufFalo-ws- 
ijmpe 



HEAVY ION PHYSICS AT RHIC 25 



0.45 




Fig. 22. a)(left) Associated charged tracks with 2 < < PTt GeV/c per trigger charged particle 
with 6 < PTt < 8 GeV/c for central Au+Au collisions at y'sjvjv = 200 GeV as a function of 
the angle AcA between the tracks in the range Inl < 1.4 compared to the data in p-p collisions 
added to the flow modulated Au+Au background 1221 b) (right) Associated charged tracks with 
1 ^ VTcu ^ 2.5 GeV/c per trigger charged particle with 2.5 < p^f 5: 4.0 GeV/c after subtraction 
of flow-modulated background. The dashed (solid) curves are the distributions that would result 
from increasing (decreasing) the flow modulation by on e unit of the systematic error; the dotted 
curve would result from decreasing by two units 11Q2I Note that only the jet correlation, after 
background subtraction, is shown. 



flow-modulated Au-|-Au background. (Note the oflFset zero.) For the p-p data there 
two peaks, a same-side peak at = where associated particles from the jet 
cluster around the trigger particle and a peak at = tt radians, from the away 
jet. For Au-|-Au central collisions, the same side peak is virtually identical to that 
in p-p collisions, while the away peak, if any, is masked by t he modulation, and, 
in any case, is much smaller than observed in p-p collisions l^SI xhe 'vanishing' of 
the away jet is consistent with jet quenching in the medium due to energy loss — the 
away parton loses energy, and perhaps stops, so that there are fewer fragments in a 
given pt range. The fact that the number of associated particles in a cone around 
the trigger particle is the same in Au-f Au as in p-p collisio ns is a strong argument 
against hadronic absorption as the cause of jet quenching l ^^^ i Since all hadrons 
would be absorbed roughly equally, the associated peak would be suppressed as 
much as the inclusive spectrum in a hadronic scenario, which is clearly not seen. 
The only escape from this conclusion is if the partons or hadrons were so strongly 
absorbed in the medium that only jets emitted from the surface were seen. Of course, 
since pions at mid-rapidity with pT > 1.4 GeV/c, jt > 10, can not be resolved 
before 14fm (— 7t/5t/™o): due to the uncertainty principle, they are formed by 
fragmentation outside the medium, even taking account of the flow velocity. Thus 
hadronic absorption in the medium is not possible for pions. 

The away jet reappears if the transverse momentum of the associated charged 
particles is lowered to the range 1 < pT^ < 2.5 GeV/c for trigger charged particles 
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with pt, in the range 2.5 < pr^ < 4.0 GeV/c (Fig^pj^^^^. In the most pe- 
ripheral coUisions, the shape of the trigger and away jets looks the same as in p-p 
collisions (Fig. [22k). However with increasing centrality, the away jet becomes much 
wider (Fig. [22b) and possibly develops a dip at Aip — tt. Since the outgoing partons 
travel much faster than the speed of sound in the medium, it has been proposed 
that a sonic-boom or mach-cone might develop, as suggested by the dip^^^^. There 
are many other ideas to explain the apparent dip, not the least of which is to get a 
better understanding of how exactly to extract the flow effect. 

Study of jet correlations in A+A collisions is much more complicated than the 
same subject in p-p collisions and one can expect a long learning curve. The next 
step in these studies is to measure jet suppression and correlations as a function 
of the angle (A0) to the reaction plane and centrality in an attempt to separate 
the effects of the density of the medium and the path length traversed. For a 
given centrality, variation of gives a variation of the path-length traversed for 
fixed initial conditions, while varying the centrality allows the initial conditions 
to vary. Although similar in spirit to a V2 measurement, this is different in detail 
since measurement of _Ryiyi(A(/>) is an absolute measurement, while V2 is a relative 
measurement, and a first look has already yielded interesting results!^. 

7.2. The smoking gun? 

The jet suppression observed at RHIC is unique in that it had never been seen in 
either p-l-A collisions or in A-t-A collisions at lower ^snn and it probes the color 
charge density of the medium. Many questions and unsolved problems remain which 
are under active investigation, but this effect comes closest of all, in the author's 
opinion, to meeting the criteria for declaring the medium a Quark Gluon Plasma: 

• There is no such effect in p-|-A collisions at any y'SNN 

• It is not the 'ordinary physics' of A-l-A coUisons since it only occurs for 
ys]^ > 30 GeV 

• In all discussions of the effect, the operative 'charge' is color and the oper- 
ative degrees of freedom are quarks and gluons. 

8. J/'S' Suppression 

Although J/^f suppression was considered to be the 'gold-plated' signature of the 

QGP and served as the inspiration for much theoretical and experimental work, 

the interpretation of the J/vP suppression observed at the GERN/SPS "^"^ in 

terms of a QGP w as controversial, and th e effe ct at the SPS could be explained by 
hadronicll™!] 

and even thermal models 1^^. The clincher would be the measure- 
ment of the J/^* at RHIG. All the initially produced J/^*, the ones not suppressed 
at GERN, would be totally suppressed in the much hotter denser QGP at RHIC, 
which would prove that the J/^E* suppression at CERN was indeed the result of de- 
confinement. However, in the ensuing years a 'nightmare scenario' developed when 
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Fig. 23. J/^ suppression rel ative to p -p collisions {Raa) as a f unction of centrality (Afpart). a) 
(left) Measurements at RHIC liMiiS] and at the CERN/SPSiM b) (right) Comparison of RHIC 
data liiS! from a) to predictions of a model 111?! with (regen.) and without (direct) recombination. 




it was real ized that if a QGP were indeed produced, the thermal c, c quarks would 
recombine | 10^ | 1Q HI to form J/'^. Thus, if the J/^ suppression were the same at 
RHIC as at CERN, this would imply that RHIC, not CERN, had discovered the 
QGP since all the initially produced J/^*, which were suppressed at CERN by what- 
ever mechanism, would be totally suppressed at RHIC, leaving only the thermal 
J/^f produced by recombination in a QGP. The nightmare is that nobody would 
believe this explanation. Incredibly, this is exactly what happened (see Fig. [23| . 
the J/^ suppression, expressed as Raa turned out to be the same at RHIC as 
in the famous SPS measurement 123. ^^^^^^^^ 

This effect is illustrated quantitatively in Fig. 123b . Modelsl ^" l ^"4n(),m 2^hich 
reproduce the SPS J/^* suppression with or without a QGP predict a near total 
absence of J/^E* at RHIC beyond 150 participants without recombination. With 
reco mbin ation turned on, the RHIC data are reproduced, with one notable excep- 
tion -^^hjch predicts a larger Raa for J/vf at RHIC than at the SPS. In fact, 
a J/\l/ enhancement would have been the smoking gun for the QGP. It will be 
interesting to see whether this occurs at the LHC. 

Does the agreement of the CERN/SPS and RHIC data for Raa in J/* pro- 
duction eliminate J /'^ suppression as a signature of deconfinement? Is it possible 
that the J/4' (c, c) is no different in its QGP sensitivity t han the 4> (s,s)-meson? 
Rec ent i ncreases in the predicted dissociation temperatures ^1^, give a possible way 
out ESI. gatz has proposed that the Xc and the ^E* were suppressed both at the 
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SPS and at RHIC, but in neither place was the direct J/^E* suppressed, since a tem- 
perature 2Tc was not reached. However temperature sufficient to melt the J/^ 
should be reached at the LHC. Does this mean that we have to wait until 2009 to 
prove or disprove an idea proposed in 1986? Fortunately there are other tests to be 
made on the RHIC data. If the Xc or were to be observed at RHIC, that would 
settle the issue. For recombination to be true, J/^I^ flow should be observed 
and both the rapidity and the pr distributions of J/ \l/ sh ould be much narrower 
due to recombination than for directly produced J/^* i 

One thing is perfectly clear from this discussion: the claim of the QGP 
discovery from J/'i suppression at the CERN/SPS was, at best, premature. 

9. Conclusions 

The medium produced at RHIC is not the expected gaseous Quark Gluon Plasma. 
It is a strongly interacting liquid of bare color charges, perhaps a perfect fluid. 
The medium flows and strongly absorbs 'colored' objects, including embedded light 
and heavy quark probes from hard-scattering, quite unlike anything observed at 
lower cm. energies. However the reaction of the medium to the J/^* seems to be 
the same, in fact nearly identical, to the effect observed at lower ^/snn- There are 
also many other mysteries to be explained such as the baryon anomaly and the 
wide and apparently split jet-correlations, possibly indicative of a Mach-cone-like 
reaction of the medium. Nevertheless, it is clear that we have moved from period 
of discovery to a period of characterizing the properties of the medium which to 
quote my colleague Ed O'Brien, "might exceed the proton in its importance to the 
understanding of non-perturbative QCD" . 
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